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MULTIVARIATE APPROXIMATION METHODS AND 
APPLICATIONS IN GEOPHYSICS AND GEODESY 


Marle-Jeanne Munteanu 
PREFACE 




This is the first report in a series which is Intended to be written by the 
author with the purpose of treating a class of approximation methods of func- 
tions in one and several variables and ways of applying them to Geophysics and 
Geodesy. 

The first report is divided in three parts and is devoted to the presentation 
of the mathematical theory and formulas. Some of my unpublished results are 
being summarized, together with my recent communications and well known re- 
sults of other authors In the field. 

We try to give a unified feeling about the topics presented not always insist- 
ing about certain details which can be found easily in the other author's papers. 

We discuss various optimal ways of representing functions in one and sev- 
eral variables and the associated error when we have information about the 
function such as satellite data of different Idnds. 

The framework chosen is Hilbert spaces, which are very important for 
their simplicity in obtaining optimal solutions, for their links with the statistical 
interpretations and inverse theory, 

Experirpents have been performed on satellite altimeter data and on satellite to satellite 
tracking data. In both cases the results are very satisfactory and they are going to be 
presented in the next report. 
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MOLTIVARIATE APPROXIMATION METHODS AND 
APPLICATIONS IN GEOPHYSICS AND GEODESY 

Marie-Jeanne Munteanu 

INTRODUCTION 

The purpose of the present survey is to make the reader aware of the existence 
of the various representations of a function of sevex’al variables and a common 
source of their derivation. 

I will discuss the numerical experience obtained with them in the satellite 
data analysis in a second report. 

The main theme of the class of methods is the representation of a function 
in a Hilbert space as a solution of a boundary value problem for a differential 
equation using Green’s function. One part of the solution can be viewed as an ap- 
proximation of the function and the second part as the corresponding remainder. 

The choice of the differential equation is very important and should be re- 
lated to the specific application. 

In the next I’eport we will discuss such a choice in connection with the sta- 
tistical interpretation and with the linlcs between different methods such as col- 
location, kernel functions, spline functions and finite elements. 

There is a wealth of connections between different approaches which should 
give a perspective on the choice of the appropriate method for the specific problem 
and also to lead us to the desired interpretations. 

The importance of the representations derived in this report is that: 

1. They furnish an approximation of a given function and the corresponding 
remainder. 
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2. An Inner product for the construction of a Milbcrt space witli reproducing 
kernel is provided. 

3. Hilbert spaces witli reproducing kernels generated in this manner are very 
interesting in optimization theory, namely numerically efficient optimal solutions 
are obtained via projection on liiffcrcnt subspaccs. 

4. Computationally very efficient optimal numerical formulas for interpola- 
tion, filtering, differentiation, integration etc. together witli their error representa- 
tions are obtained very elegantly. 

5. Error estimates for the optimal formulas in Hilbert spaces arc given by 
tlic hypereycle inequality. 

6. The reproducing kernel of a Hilbert space provides us with a covariance 
function an?j therefore leads us to statistical interpretation. 

7. Hilbert space treatment seems to be an alternative to the statistical 
approach. 

8. Generalized inverse problems wJiich are so vitally important in geophysics 
have many times as a framework the Hilbert spaces approaches. (Sec Bachus and 
Gilbert papers.) All these ideas will be discussed in detail in future reports. 

Concerning tlie numerical c.\j)eriencc with these methods, together with Carl 
Sicfringl have experimented with the tensor product representation and more 
specifically with the tensor product of univariate B-spline functions. 

In the application of the tensor products one needs a first stage, namely gen- 
erating the values of the unknown function on tiic rectangular grid with, forex- 
ample, least squares and orthogonal polynomials. 

We are applying this scheme to sea surface topography. 


Qn the test data the tensor product of splines works outstandingly well, the 
error associated with the method is extremely small. What is extremely inter- 
esting about the tensor product of B-splines is that the univariate B-splincs arc 
functions with support of small size am/ outstanding numerical properties so that 
the above two dimensional sclieme can be viewed as a v;ry efficient finite element 
formula. This finite elemefit scheme can be used not only in the best approxima- 
tion context such as optimal filtering or interpolation but also in optimally solving 
partial differential equations as well as optimization problems of ail kinds. In geo- 
physical problems such as for example modelling the SAN ANDREAS fault or in the 
variational approach of whole-earth dynamics, etc. we can use these functions as 
trial functions. Our numerical experimentation witli satellite data such as altimeter 
data is going to be included in the next report and hopefully be presented at the AGU 
meeting this spring. Of cour.se the nature of for example altimeter data also calls for 
tlie application of other schemes such as blending .spline functions. Tins is the next 
step in sequence. The reason is that tensor product rpline functions interpolate data 
on a rectangular grid, while blending spline functions interpolate data along curves 
in space which in our case arc the satellite tracks. 

In addition Jf wc use B-spIine and least squares in the blending schemes one can 
compress data significantly. We have experimented very succesfully witli satellite to 
satellite tracking data, compressing around 300 data to 32 coefficients. 

'file FORTRAN subroutines used and the graphs obtained are available and 
will be included in a note together with all the other experiments. 
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PART I 


The first part is concerned with the presentation of a general method in order 
to obtain representations of functions in several variables. This uses, essentially, 
the form of the solution of a boundary value problem for an oi'dlnary differential 
equation and the isomorphism between the lattice of a collection of linear projec- 
tion operators and the lattice of corresponding image spaces. New representa- 
tions, well known formulas such as the Taylorian representation of D.D, Stancu 
[29] and Sard [2o] as well as the vei’y interesting and elegant results concerning 
multivariate functions obtained by Gordon [19] , Gordon and Birkhoff [7] are ob- 
tained using intuitive diagrams. The reproducing kernel for a Hilbert space of 
bivariate functions which have Ta 3 'lor representations can then be constructed 
(see Nielson [24] , Mansfield [21] , and Ritter [25]). 

One can generalize their method for more general cases than bivariate spaces 
of functions possesing Taylor expansions and their variations. Once the reproduc- 
ing kernels are at hand the question of optimal approximation of linear functionals 
in Hilbert spaces is reduced to requiring that the approximation be exact for the 
representers of approximating functionals. 
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PART II 


In the second part approximation 1‘ormuias and tlie corresponding remainders 
are characterized as solutions of certain partial differential problems. These 
formulas furnish representations of functions belonging to o. certain space of bi- 
variate functions and provide in some cases an inner product. Tiie space c.m be 
completed in the sense of this inner product, reproducing kernels kernels can be 
constructed tind consequently optimal approximations of linear functionals and 
other very useful quantities in numerical analysis can be derived. One can de- 
idve an entire variety of differential problems giving intermediate approxima- 
tions and the corresponding remainders, maximal and minimal approximations 
being the two extreme cases. 

The isomorphism between the lattice of a class of differential problems and 
the corresponding lattice of linear projection operators which are solutions of 
the differ’ential problems is used in the main. If, in addition, the isomorphism 
to the corresponding lattice generated by the projectors is used the theory of 
approximation of multivariate functions can be unified in a beautiful manner. 

At the end of the second part of the report the tensor product splines of 
Ahlberg, Nilson and Walsh, ^Jordon's bleading splines interpolating a network of 
curves and generally splines in several variables interpolating hypersurTaces 
and corresponding remainders are introduced. It is well known that blending 
splines in two variables have been used with great success in problems of 
computer-aided design, airplane fuselages, automobile exteriors, etc. Gener- 
ally speaking these multivariate splines cmi be used in the problem of modeling 
smooth surfaces. 
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Tiie major role in the construction of these multivariate spline functions is 
played by univariate spline functions. Using Gordon's language the univariate 
spline functions generate a distributive lattice of projectors which has a "maxi- 
mal" projector and a "minimal" projetitor. For the bivariate case the "minimal" 
projector is the tensor product spline and the "maximal" projector is the blend- 
ing spline mentioned above. 

The lattice of corresp<^nd:rg remainders gives the errors associated with 
the maximal, minimal and generally all intermediate approximations. 

The above mentioned isomorphism between lattices suggests the Euler- 
Lagrange differential equations satisfied by these splines. Consequently, the 
minimization properties of these functions can be derived. 

In a later report several approaches to defining multivariate splines will be 
presented. Of course, one important approach will be the variational one, start- 
ing with the minimization property of spline functions, deducing the Euler -Lagrange 
differential equations and finally obtaining the analytical fr-'m of certain types of 
multivariate splines. 

1. REPRESENTATIONS OF FUNCTIONS JI’l SEVERAL VARIABLES 

1.1 The Bivariate Case 

1.1.1 W. J. Gordon's Maximal and Minimal Approximations 

Consider the space C»’[I'] , the set of real-valued functions, defined and con- 
tinuous on the interval I' »= [0,1],^ such that all derivatives of order less than or 
equal to p are continuous o.c I'. Similarly we define the space C [I"j where 
I" = [0,ljy. Consider also the space C’'' ER^J of bivariate functions, defined 
on R^= I'x I", real -’Valued, and having derivatives 
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(i.lj 

y)fC’’(R^j ) OSiSp, O^jicj, where 

pO.i), 

Sx'Byi 

Let D^ and Dy be Ihe lineai' differential operators 

= Hp (x) ~1 + . . . + n j (X) n^(x) , 

dx** 


“v= 


• +^i(y)^ + t>o(y). 


defined on C” [I*] , C*' [!"] respectively. We assume aj(x), (y) are x‘eal and 
continuous on their domains of definition.s and a (x) ^ 0, b (y) 0. 

Denote by and Ny the respective null spaces of the operators and Dy, 
Tliese are obviously subspaces of dimension p and q respectively. 


Let , 1* 0,p-l iind nij', j *= 0,q~l be continuous linear functionals on 
C’LI’], respectively, which form the bases for the dual null spaces. The 


functionals i= 0,p-l; lUj, j * 0,q~l are of the following form; 


^ ^ 5>*“1 

Au.f<^t)+E f -flkWfSdx.fsC'U'). 

k»0 i‘=° 


Vi ^ (fc') r (k') 

in'(:E) = 2] L, 6(x^) + 2] ak'(y)B(y)dy, geC'’ [I"]. 

k'-O t'»0 k'nO •'l'* 

where Aj^^, are constants, the functions ^^{x) and are piecewise 

continuous on I' and I" respectively, and y^, ax*e points in I' and I" respectively. 

Denote by -f j, i = 0,p-l and , j = 0,q-l the respective restrictions of the 
functionals -t'; , i « Oip-l and m I , j = O';^ to the subspaces and Ny . We know 
that there exists a unique base (x), i = 0,p-l for Nx and a unique base \fj. (y), 

3 * O^q-i for Ny sucli that 
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I ~ ^ik' k = 0 , p 1 , 

'"jC'^'k) = S|,i. k, j = 0,q-l. 

Lamina 1; Every function f(x)£C’’ [!'] can be represented in tlie following form: 


p-i 




1.0 


Gj(x, t)Dj(f)dt 


( 1 ) 


where G (x,t) is the Green's function of the differential problem 


D^f = c 


^ p /rx n • r where feC'"[l'] — E s_^g eL^ tl'] . 

(f) = 0 1 = 0, p -1 

(see [4] , [10]). Similarly all g(y)e C’’ [!"] can be written as 

<1-1 

g(y) = V mj (g)i//. (y) + G 2 (y-‘')D„ (B)du 

j.o •'i" 

G_(y,u) is the Gx'een's function of the differential problem 


( 2 ) 


DyE = b 


gecni' ] 5-—heL^ [I"3 


,n>j <e) = 0 j = 0, q-1 

Define the "partial" operators Lj, i = 0,p-l and , j = 0,q-l on the space 
C’’’‘’[R 3 such that 

L.(F) =.ej[F(x,y*)] for ally* I", 


F(x,y) C’’-'' [R^] , 


Mj (■’ . = nij [F(x*,y)] for all x* I'. 


Denote by D and D 




the extensions of the operators D^^ and Dy to the space c'’" '' [R„]. 
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Theox^em i Eveiy F(x,y)eC‘’*‘' [R^^] can be written in the form 

P-i - 

F(x,y)= ^ L,(F)0j(x)+ G,(x.t)DjF(t,y)]dt=:Pj(F) + Rj(F) (3) 

, iaO *'l' 

and 

pj r 

F(x.y)= 2 ] M^(F).Aj(y) + G2Cy,u)DjF(x,ii)]du = P2(F) + R2(F) ( 4 ) 

J=o ■'i" 

(see Gordon [17], p. 17). 

It is obvious that all functions F(x,y) eN HN can be written as 

X y 

p-l q-l 

F(«*y)=E I]LiMjCF)0i(x)./-j(y)=PjP2(F). (5) 

i eO J bQ 

Using the results given in [17] , [7] , [19] were able to assert that all P(x,y) 
c N'flN' can be represented as 

X y 


t(x,y) = 


G, (x,t) G^Cy.u) Dj DyF(t,u) dt du 


(6) 


where and are the complementary spaces of and Ny , respectively. 
Suppose now we wish to obtain our function in the form 

F = P(F) + R(F) 

whore P(F) is the component of the function belonging to the subspace N,nN„ and 

X y 

satisfying the Interpolation conditions: 

Lj(F) = Li(P), i = 0,p-l, 


Mj(F).= Mj(P), j = 0,q-l, 

and R(F) if the corresponding error, i.e. the component of the function belonging 
to the complementary subspace of N UN . 

Let us first make two remarks about the n'^'’ order. These notes will be 
used constantly in this chapter. 
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1. UBinK Dif'udonne's theorem (in [13] ,p. 123) we can obtain F as the sum 
of two coiii|)onent8, one in the subspace and the other the corresponding error. 

2. There exists an isomorphism between the lattice of the subspaces and 
the lattice of the corresponding representations. (See Gordon il9]). 

Theorem 1: Kvery function F(x,y)£ C’*'*’ [R^,] admits the following development: 

F(x.y) = P(F) f R(F), 

where 

p-I q-l p-l q-l 

P(F)= 2] L^(F)<?>.(x) + ^ M,(F).//^(y) - 2] 2] (V) 

•■0 )m0 iaO JaO 

= P,(F) + Pj(F)- P,P,(F). 

*^(F)= jj i,(x.t)Gj(y.ii) D,D„F(t.u)dt du. 

^11 

By using the following diagram, remarks 1 and 2, and formulas (3), (4), (5) a 
and (6), we arrive at formula (7). 

We can e.isily verify that P(F) satisfies all the necessary interpolation 
conditions. 

N.n = (N, uNy)' 


This diagram suggests that we compose Pj (F) and PjfF) to obtain P(F). 
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Denote this opcvnlion by u) , then 

P(K) = (P,® P2)(F) = P,(F) + P3(F)-PiP^(F). 

In the terminology intvodncod by Gordon ([19] p. 234) P{F) is called tlio miLxinial 
approximation. 

We now look for F in the form 

P(x,y) = p(F)+ r(F). 

whore p(F) is tlie component of F in the subspace NON satisfying the inter- 
potation conditions 

LiMj(F) = LiMj [p(F)] , i O.p-l, j = OTq^. 
and p(F) is tlio corresponding error, i.o, the component of F in tho complementary 
subspace of N . fl N. 


Theorem 2 ; Any F(x,y)£C'*’‘’ [Rj can be written in tho form 

F= p(F) + r(F). 


where 


p-i n-i 
IkO i-O 


f(F)^ j G,(x.t)D^F(t,y)dt .|. f GjCy.u) D^,F(,u) du ... 

I ' A " 

G,(x.t)G 2 (y.u)Dj D/(t,u)dUlu = 

z:R,(F).t.R2(F)-R,R2(F)= CR,(dG2)(F). (8) 

An immediate justification of the formula is suggested by the following 
diagram: 
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generiliization of a development of 


Remarks: 

1. Theorem 2 can be considered as a 
Stancu(see [29]), 

In fact, if one particularizes as follows, 

r = [0, !],!'= [0, 1] . 
L,M.(F) = D^Vo.i)). i = on^^, 
one obtains the taylorian development of Stancu: 


P-l q-l 

X* yJ 

- y --'0,0) + 


I’-Cx.y) = T' TT ^^Co/o)+ L_ r' . p-i (P.O) 

uo J' (P-1)! J t); ‘D^(t.y)dt 

J 'y-^)rD^‘’*'’V(x,u)du + 


r'r' 

4 


(y-o)2 


q-l 


(p-1)! (g-1)! du 

Obviously p,., satisfies tbesecessary interpolation conditions. Condon 
calls p(f) the minimal approximation. 

The diagram in this case gives us the errors K.<., end E,,P, „f the opera- 
tion of composition, then r(F), 

We can obtain another expression i ■ 

expression foi r{F) by making another diagram. 
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Theorem 3 ; Every function F(x,y) t C’*'*’ (II^) can be written in the following 
form: 

F(x,y) r P,P,(F) ♦ P,R,(F) + PjR,(F) + R,R,(F) r 


p-i <j-i p-i - 

f52(y.‘OD„ lL^(F)]du 

• ■0 J >0 i.O 

+ >_]Vi(y) f C,(x.t)D, lM,(F)]<it ♦ ff G,(x.t)Gj(y.ii)D,D„(F)dtdu 

i-o *r ••^R 


( 10 ) 


We justify this formula in a manner ick>ntical to the preceding one but we 
begin with the following diagram: 



N' n N 
» y 


n.on; 


For example, for the shaded part we have the following conditions; F(x,y) 
£ N and F(x,y) e N' . Then we have the representations 


p-i 


F(x.y) 


F(x.y) = 2] L,(F)(^^(x) 

i.O 

■I 


Gj(y,u)D^, F(x,n)dii 


(H) 


(12) 


By substituting (11) into (12) we see that all functions F(x,y) t N n N' can 

* y 

Ije written in the form 


F(x,y) 


P-« . 

= 21 0.(x) J o, 


(ydi)D„ Ibj (F)l dll. 
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Similarly all F(x,y) e N n N' c{in bo writte.u as 


1-1 . 

F(x,y) 2] \^j(y) G,(x,t)D, [Mj(F)l 

j»o ‘^r’ 


dt. 


In particular sines we have already proposed Theorem 2 we obtain a taylor- 
itm development whose expression Is different from the rest. 

Formulas (7) and (LO) were deduced by Gordon in a different way (see [18] 
p. 939, [17] p. 23, 2.5). 

We have demonstrated a method in the very simple case of bivariate func- 
tions and presented for them the maximal and minimal approximations, Tliis 
discussion is an introduction to the case of several variables and to tlie subse- 
quently more complicated representations. 

The above formula gives us a very interesting way of introducing a scalar 
product for the bivariate function spaces. For a particular simple choice of the 
operators and the functionals we obtain one of the Sard spaces [2C] which is very 
often used in numerical analysis. 


1.1.2 Geiieralizatioas of the Development of A. Saxxl 

Let = I' X I", (a,b) e and p, q are non-negative integers such that 
p+q=n. A, Sard in [26] introduces the space (a,b) * ^ (R^; a,b) 

as the collection of real-valued functions f(x,y) defined on R^^ and having derivatives 

(p.q) 

Bf(x.y) (x,y) £ 

Df(x,b) xe I' j < q, 

(i .n-J ) 

D^a-y) yel" i <P. 

which are continuous. Assume p + q * 2m, a» 0, b * 0. 
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Theorom 4 ; Every function f(x,y) « ^ can be written in the form: 

f = A i U « C i D (13) 

where 


A = 



t *|<2a 



(y-u) 


U.lm-l) 

Dt (0,u)du, 


!<■ •'O 




( 2m-j>l )! 


;D^(t‘Wit 




(x-oT'* (x-or* (y-u)r* 

(m-l)! (m-1)! (m-1)! 


(«.■) 

Dtt.u) 


u)dt(iu. 


To begin we will write out in detail the first sum. 


where 


■-I 


w-l 


imO 

■-I «-l 


*■'1 Ifrfr 


1-0 


iaO JaO 


■ *l 

m-l 

2iK-l-i 


= L 

z: 

iaO 

imO 

i-« 

■ -1 

m -1 

2«-l-l 

I V, 


I 

1-0 

1-0 

i mm 


* vl (,*• I ) 


Construct the following diagram: 


(KerD;)' (Kerl)/*)' 
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Denote by 


Di tr — , d‘ =ii, i = 0.2m, 

"ax' ay' 

Also denovC by S[x] , Sly*”, . . . , y^'fl-i] , the subspaces generated by x and y™, 

. . . , y2>n-i j respectively. 

It is clear that all f(x,y) t Ker D™ n Ker D™ admits the representation Aj, 
all f(x,y) e Sty'] fl Sfx'^^’, .... x^'"“^”'] , j « 0,m-l the x'epresentation A^^ 

J = 0,m-l and all functions f(x,y) e sLx'l n S[y'“, . . . , representa- 
tion A 2 , i , 1 0,m-l, 

Search the subspaces now for the one wliich has the B representation. Fol- 


lowing the same diagram we prove that all f(x,y) which ax’e elements of the sub- 
space {S[lJ n {Ker dJ"")'} U{S[x] n (Ker )'} U . . . U (slx'"] (Ker 
can be written in the B form. 

Actually for all f(x,y) e S[x'] fl (Ker Dy"'“' )', i« 0,m-l, we have that 


f(x.y) = 




1 / 2i 

(y-'-0+ 


(2m-i-l)i 


f(0,u)du, i = 0,m-l. 


Similai'l 3 ' we easily see the C representation corpesponds to the subspace 
{S[l] n (Ker dJ™)'} U {S[yl n (Ker U . .. . U (Sfy"*] n (Ker D™^^)'}. 

Finally we Icnow that all f(x,y) e (Ker D'")' fl (Ker Dy)' admit the D repre- 


sentation and fox'mula (13) is thus proved. 


Remarks 

1. Formula (13) of Theorem 4 constitutes the taylorian developzxxent of Sard. 
The author has obtained this representation in u completely different manner. 
(See [2G]). 
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2. The method, the one by which we have just constructed the formula of 
Sard, permits us to obtain more general representations possessing similar 
characteristics. We Imve presented the method in the trivial case. We are able 
to obtain in an analogous manner more complex representations of the form 
(notation is obvious) 

f(x,y)= + XI f G2m-i Cy.u)D^(t!u) fL,(f)]du + 

i»J<2ni Km “'0 

j<m *'0 



(14) 


We have replaced the taylorian functionals by the continuous linear functionals 
LjMj , I + j < 2m. 

This representation Is considered as a generalization of the taylorian de- 
velopment of Sard. 
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2. Decomposition o£ Linear Differential Problems 

In the second p*irt of the I’oport we will obtain representations of multi- 
variable functions starting with appropriate decomposition of certain linear dif- 
ferential problems. We also find again certain developments obtained in the 
preceding part of the report. We organize the miaterial in tlie following way. 

Two differential problems, for the bivariate case, are treated in the first 
section and we obtain the maximal and minimal appro-'^imations and the corre- 
sponding ei*rox’S. One ease in particular permits us to concretely illustrate a 
result given by G. Birldioff and W. J. Gordon [7] , 

In the second section we study two differential problems for the four variable 
case whose solutions can bo interpreted as defining an interpolating function for 
a system of given curves and for a system of given surfaces, respectively. We 
also study the corresponding error’s. As one last example, for the four variable 
case, we intend to illustrate a partial symmetric intermediate approximation. 
Finally we consider the general case of n variables as providing the inter- 
polating function of certain hypersurfaces in j variables, j = T,n-i and the cor- 
responding errors as solutions of certain differential problems. For the n 

variable case we are able to obtain all tlie total S 3 annietric intermediate approx- 
imations and the corresponding eri'ors. 

From the point of view of the theoi’i' of differential equations we obtained 
tlio solution to a certain class of differential problems by systematically using 
Ihe same method based on the appropriate diagrams, suggesting the mentioned 
isomorphism between corresponding distributive lattices. 
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2.1 Bivariate Case 


Wo adopt the same framework as found in part one. 
IjCt the non-homogeneous differential problem be 
'j9^J!)yF(x,y) = e(x,y), 


Li(F)- fj(y). 
Lmj(F)=:e,(x). 


i = 0, p-1, 


j=0,q-l, (17) 

the symbolism is the same as in part one, e(x,y) e C°(R^), fj(y) e O'* [I"], 


g. (x) e C [!'] , i * 0,p-l, j 0, q-1. 

Suppose that the compatibility conditions 


LjCEj ) = Mj (fj ) . I = 0. P-1 . j = 0,q-l , 
are satisfied. We can then assume a solution of the form 

F(x,y) =:Fj(x,y)+ F^Cx.y), Fg C'’*'’ [Rj 
whore Fi(x,y) is the solution to the differential problem 

'f,i9yFi(x,y) = 0 


i = o. p-i. 


-Mj(Fj) = g 3 (x), j =0,q-l, 


and F^ (x,y) is the solution of the problem 




L. (Fj ) = 0, i = 0, p-1 , 


_M.(F 2 ) = 0, j=0,q-l 


We then have from the first part that is 


FjCx.y) = jj Gj(x,t) Gjfy.u) e(t,u) dtdu 


(18) 


(19) 


( 20 ) 


m 


( 22 ) 
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Following tlio diagram corrosponding to Theorem 1, part one, and tl. 
relation 

+ ^KcrlSJj " ^K-srlS)^ KerjSy' 

where by we moan tho operator corresponding to all F(x,y) e [RJ , 

the component in tlio subspace Kor 1!) i'l , We decompose the second problem into 
three diiferontial problems: 

= 0, 


,Li(F„)= fi(y), i = 0,p-l, 

'S,F.2 = O’ 


.M,(F, 2 ) = Bj W. i=0, q-I, 
■»«F,3 = 0, 


^L.M^(Fi 3 ):.Lj(g.)=ajj. i = 0,p-l, j=0,q-l, 

for which the solutions are: 


p-i 


f,. 5,L‘ m &)■?>,«. 


q-1 


''.j = E h '/•i (y) ' 


i-o 


^i3 = fi fi(y) 

1=0 j=0 


From ( 23 ) we conclude 


from which it follows that 


F r: F 4- F — F 

11 ^ * 12 ^13 


F = F +F ~F +■5' 

I 11 + 1 12~ 1 J ^ + -2 • 


(24) 
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This is tho result proposed by Theorem 1, part one. The uniqueness of the solu- 
tion to problem (17) is due to the hypothesis made on the operators Lj and Mj , 
(see part one) . 


Remark: 

l^or the particular case: *= pcsp.o) ^ * p( 0 r 2 p) ^ 

L.Mj (F) = i = 0, p-1, j - £\ p-1, X,, 0 or 1, 

G. Birldioff and W. J. Gordu' i obtained the same result using an induction method 
(see [7]). 


Particular Case 

We are going to consider a piu'ticular case which appears in [7] . The solu- 
tion given by the author on p, 204, formula (19) is inexact and will be corrected 
in what follows: 


F^x,y^ = 0, 0 £ x,y 5 1 , 

F(x,0)= 0, F(0,y) = 0 
■I F(x.l) = 0, F(l,y) = 0 

FyC<.0) = BoC‘)> F,,<o.y)= fo(y). 

F,(l,y)= ft(y), 

where f. (y), qj(x) e C'LO,!] , l,j » 0,1, and we have 

i.j =0,1, 
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Tlitn the unique solution In C*’’’’ [R„I is of the form 

F'C'^.y = foC’*) + “^iCy) fi(x) + ^o(^) 6o(y) +‘/’i(x)gj(y) - 

- ^o(y) F(O'O) ~ W‘^o(^)*^i(y) F*(W) - 0i(x)0o(y) 

- 4i(s.)4 i(y)F(l.l) - ‘f2(5c)0o(y)F(OiO) ~ (y)^^.^! “ 

“■^3<’^>>"o(y)F(o:i) “03(^)^i(y)FH,’ii -02(y)0o(«)F^M5'' 

. <0.1i (0.1) (0,1) 

-‘/2(y)^i(-'^')F(0.0) ~^3(y)^oC>i)FCl,O3-<^3(y)0j(x)Fa,l), 


^ where 

0j(x) = x2(l „ x), 
stjCx) r= (1 + 2 x)(1-.x)2, 

03(X)=: (3 -2x)x2, 

are the cardinal functions for the polynomial Interpolation problem in tl»e sense 
of liermite (see Davis 1,12] p. 45). 

The case where tlie four initial boundary conditions are non-zero has been 
conslderiid by Gordon in [17] p. 38. 

Let n '• ,iow consider the following differential problem; 

'■i9,F= fj(x,y). 

J fa(x.y). 

Ll.m^(F) =: inOTFT. j = o7?pl, 

where f.{j:,y) e C°[I'] x cni"] , faix.y) e CPII'] x c°[r']. The following 
condition 

^^fjCx.y) - J5yfj(x,y) (26) 
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is assumed to be satisfied, then every solution F(x,y) is necessarily a sum 

F(x,y) = Fi(x,y) + Fj(x,y) 
where Fj (x,y) is the solution of the problem 


rs.F, = 0 




S„F, = 0 


(Fj) = a.j, i = 0,p-l, j=0,q-l, 


and Fj (x,j') is the solution of the problem 

fj(x.y). 


.L.M^(F 2 ) = 0, i = 0,p-l, j=0,q-l, 

The solution Fj (x,y) is obviously given by 

p-l q-l 

Fj(x,y) = 2] 

jaO 


(27) 


(28) 


(29) 


Looking for the expression of Fj we return to the diagram corresponding to 
Theorem 2, part one. Let T be the operator considered in the problem. Then 
we deduce that 


^(KcrT)' ”^(KorJ9^)' +^(Kerl9^)' " ^<KerS„ ) 'fl (Ker J?„ ) ' 

X y * y 

by (KerT) ' we mean the complementary subspace of the null space of the opera- 
tor T. 
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Then to find the expression for Fj (x,y) we decompose the differential problem 
(28) into three differential pi’oblems; 

Ll. (Fjj) = 0. i = 0, p-1 , 

= f2C’^-y) 


M.(F 22 ) = 0, j = 0,q-l 

= jgjjCx.y) 


l.(F 23 )-o, i = o,p-i 


Ulj<F23) = 0> j=0.q-l 

For the reasons given in part one > we then have 


n 


F, 
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F, 


[ Gi(y. 
Jq 

■=r 

“'o 


Gj(x,t) fj(t,y)dt, 
u) f2(x,u)du, 


Gj(x,t)G2(y,u) f2(t,u)dtdu. 

I ' ing calculated the relation in (30) we obtain 

F2=F2i+F22-F23. (31) 

Consider the same differential problem and work this time for the expres- 
sion of Fo from the diagram corresponding to Theorem 3, part one. Thus 

^=■2 = ^2 + + Ff <32) 


where 

(1) F^' satisfies the differential problem 


U.(f;) = o. i = o.p-i. 
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and con'OspoiulB to Uio mill spac»i< of tho oponUor 

J«0 


(2) BatlafloB tho problem 




0, j=0,q-l , 
luul oorrospontls to tho null spaco of tho oporator JDy , 

l»Q 

(y) lilvo above, la tho solatium to tho problem 


< Li(i-'3’’ = o, i = ori^, 


It Is a clear result that 


q-l r 1 

K - E •/'] (y) G , (X, t) M J [ I'l (t ,y)] cl t , 
J mO •'O 

J>-1 ,• I 

^^2 “ Z] C^aCiM') L, [fjCx.u)] du, 

I -n •'n 


UQ 
I ri 


J f Gi(x, 

*'0 *'0 


(x,t) G 2 (y,u) 19^ f 2 (t,u)dtdu. 


This Is tho result given by Thcorom 3, part ono. 

In an analogous mannoi- wo can treat tho case of throe and inoro than throe 
variables. 
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2.2 Spline FvmcUons in Sovevai Vnrinblos 

This Hoction is inoiini to bi'leOy present in a uniriod manner certain types 
of spline fiinetions in several variables. The major vole in the constinction ol' 
these spline functions is played by univariate spline functiona. 

The reasons behind plaeinp; this section at the end of this report are; 

1. Wo use tljo methods of part one, sections one and two in order to derive 
the analytic expression of the Spline funetion and the corresponding remnindov. 

2. The differential problems studied In part two of Iho report suggest to us 
the Eulor-Lagr;mge equations of tlio minimization properties satisfied by tlie 
spline funetions under consideration. 
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2.2.1 I^dcollootions 


Let 77' : 0 »= jCg < Xi < . . . < Xj,, - 1, 0 « < yj < . . . < y^, - 1 ftnd 

77"': 0*Zj<z,<...<Z|.,»=lbo three ptirtltlons of the intervals I' * C0,l] ^ • 
I"* (Oa]^., and I"'« [0,1],. 

Let E^i [V] be the collection of functions f defined on 1' which are of class 
C^’’ inside each subinterval [Xj , Xj^j] of 1* luid which have simple jump discon- 
tinuities in the (2m-2)'^'' derivative aci'oss the joints Xj , 1 < i < p - 1 . Similarly 
define the space E=';[l"] and [I"»] . 


Define '’«[!' i"] as the space of functions f{x,y) defined on I’ x r* and 

whose normal derivatives satisfy the conditions 

(1“) [I"], Vxx-e r'. 0 <M<2m-2, 

(2“) f(x,y*)e E^“; Cl'] , -Vy* e I", 0 < v < 2m - 2, 

For the trivariate case we introduce, in an analogous fashion, the space 

^ ^ partition 77' x 77"x?7"’ of K,, yields 

purallelpipeds denoted R , i = 0^*, j ® 0^, k « 0^, 

I.et the partial operators be 




Lj (f) = 'Cj [f(x,y*, 2 *)] , for all y* 2 >fcj,"x I"’, 

M. (f) = m. [ f(x*,y, 2 '»‘)] , for ni I x* 2 * u I' x I"', 

(f) = n,^ C f(x 2 .)] . f or n 1 1 x*,y* e I' X I' , 


i = 1, P, j = 1, Q, Ic = 1, R. 


'rhe functionals 'll . , 1 < i ;< P are linear, continuous, defined on E^*,’ [!'] , Jind 

>. jff 

of the following form: 

ffCxg). *£ = 0, p' 

Vf(x)€E;‘;[I'],.C,Cf(x)] ^ 

't = 0,p'. a=l,p-l. 

Similarly on tiie spaces E^|) [!"] , E*|i [I*”] we define the functionals m, . j « LQ 



Wc introduce ns In sections 1,2 and part one differential operators S , 
of degi'eo p, whore P, Q, II > p. 

Denote by the collection of spline functions related to the differential 
operator imd the functionals 1 i i < p. The dimension of is P. Let 
i the cardinal base of S, . 

A spline function a(x) t for the case of the differential operator D’’, has 
the following properties 

(i) s(x) is a polynojnial of degree less tlmn or equal to 2p-l in eacli interval 

oIp'Q, 

(ii) the first 2p-2 derivatives of sfx) lU'e continuous at the joints x, , i * Ojp^l. 
Starting Avith this unified definition wo have for the trivariate case, we can 

generate the spline functions of Gordon luid Ajilberg-Nilson- Walsh. 

2.2.2 W. J. Gordon* s Spline Punotions 

Suppose we wish to obtain an interpolation formula 

f(x,y, 2 ) S(x,y,z) + r(x,y, 2 ), (49) 

such that 

S(x,y.2)r. 0 (50) 

in the interior of the paralielpipeds i >= 0^ , j * oTT’, k » 0^, and that 
the following interpolation conditions are satisfied: 

L.tf)^L^(S), I-TTp, 

M^(0-M.(S), j«LQ. 

N^(f) N^(S), ks^ITR. (51) 

Wg denote by the adjoint form of tiic operator D. 

The spline function s(x,y,z) was defined in a different way by Gordon [ISj. 

To obtain tlie analytic expression of s(x,y,z) and v(x>y,z) wo construct the 
following diagram: 
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Via iinalogouB reasons as those of sections 1,2 and par* 1 we can write the 
iollowing formulas: 

p 0 R 

S(x.y.z) = + I]N^(f) X^(z) - 

IbI 

P 0 PR 

‘11 "L 2]t.N,(f).;.,(x)X,(z) - (52) 

iatjal ialksl 


OR r 0 R 

)b 1 k>t IbI l.t kBi 

r = jjj Ci"(x. t)Gj(y, 11)63(2, v)f)* 4 ), f(t, II. v)dt(' idv (53) 

where Gj(x,t) is the function which appears in the error formula for a univariate 


spline 


f (x) r S(x) + r(x) 


S(x) = 2^ (x) 


r‘ f 

(x) = I |G,(x,t)- 2 ] (G(x,t)]j f(t)dt = 

•'0 IbI * 

= [* C“(x,t)D*Dj f (t)dt . 

•'ft 
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Gj(x,t) is the corresponding Green's function of tiie boundar}' value problem 


■DXf(;0 = o(x) 

(a) 

f(0) = 0. I 

— — 

,f(l) = 0, a = o.p-l. 


f(x) € C*’’ [I 


D*D 

X 


In tile sum whicli appears in the integral, we consider all the functionals , 
i = l,P,with the exception of those in the boundary value problem (57). (See 
Birldioff and de Boor [6]). 

The functions G 2 (y,u) and G 3 (z,u) are similarly defined. 

The analytic expressions of tiie spline function s(x,y, 2 ) and the error r(x,y,z) 
were obtained in a different way than Gordon's, in [18] , for the bivariate case. 


Minimization Propert}'^ 

Gordon's spline function minimizes the pseudo-norm 

f(x.y,z)]‘dxdydz, (58) 

with respect to all functions wiiich are in the space [R^i] , and which satisfy 

tiie interpolation conditions (51). (See [IS] , for the example). 



2.2.3 Multi-vainable Spline Functions Introduced by 
Ahlberg, Nilson and Walsh 

We are looking for an interpolation formula 

f(x,y,z) = Sj(x,y,z) + i(x,y,z), (59) 

such that Sj(x,y,z) satisfies the equations 

ft SjCx.y.z) = 0, 

fty Si(x,y,z) = 0, 

ft S,(x,y,z) = 0, (60) 
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in each of the parallelplpeds Rj j |j» with the follo\vdng interpolation conditions 

L^M.N,,(Si) = L,M^N^(f). = i = ITp , j = i 7Q. k=l7R (61) 

This definition is given in [1], The authors have not studied the corresponding 
error (see also [9]). The diagram resembes that of part one section two but we 
replace Nj^, Ny, Nj, by Sy, Sj.. We have the formulas: 

p <3 R 

S,(x,y,z) = ^ £ E LiMjN,(f)0i(x)0j(y)X,(z). (62) 

Ul jsl ksl 

p <3 f 

r(x,y,z)=2^ 2] (j^)'/'j(y) G3(z,v)J9^fi^ [LjMj(f)Jdv + 

ial j = i ‘^1'" 

*L Z;^j<y>V‘=) f GU’‘.‘)«r«, [Mj\(f)jdt + f; f]0.(x)x,(z). 

j = l k=l ‘=1 kal 

(* G°(y,u)J!)X[LiN,(f)]du-2]0i(x) [f G“(y,u) G 3 (z,v)igXfivX fL.(f)] dudv 
- 2]^ >/'j(y) ff Gj(x,t)G3(z,v)J£>^J5)^rt^[M.(f)]dtdv ~ 

“ Z] (*f Gi(x,t)G°(y,u)i!)*J!)^J!)^f^^^ [N^(f)] dtdu 

k=ri 

G°(x,t)G°(y,u)G3(2,v)i)fj9^i9X^vX If(t,u.v)J dtdudv. (63) 

Minimization Property 

The spline function Sj(x,y,z) minimizes the pseudo-norm of (58) with respect 
to all the functions in the space which satisfy the inlterpolation con- 

ditions (61). (See [!])• 

REMARKS 

We can similarly introduce spline functions in n variables interpolating 
hypersurfaces in j variables j * l,n-l. Using univariate spline functions we can 
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generate a lattice of projectors whose "minimal" projector will be the tensor 
product splines and whose "maximal" projector will be a spline function in n 
variables, interpolating a network of hype rsurf aces in n - 1 variables. Of course 
of great interest will be all the spectrum of intermediate splines interpolating 
hypex’surfaces in any number of variables from one to n - 1. It is unnecessary 
to mention the great importance to have at hand tlie lattice of corresponding re- 
mainders for numerous problems In numerical analysis. 
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